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Abstract. - I construct a well-defined expansion in e = 2 — d for diffusion processes on small- 
world networks. The technique permits one to calculate the average over disorder of moments 
of the Green's function, and is used to calculate the average Green's function and fluctuations 
to first non-leading order in e, giving results which agree with numerics. This technique is also 
applicable to other problems of diffusion in random media. 



The small- world network [1] has served as a fundamental model in the field of networks [2]. 
However, the problem of averaging over the possible different random connections in the small- 
world network is severe in low dimensions: a study [3] of the properties of even the simple 
problem of diffusion on the one-dimensional network leads to a difficult problem that, thus 
far, has only been tackled approximately. 

The physical reason for this problem is a breakdown of mean-field theory [4] in dimensions 
d less than two, and the emergence of strong site-to-site fluctuations of the Green's function, 
so that the properties of the system cannot be represented by simply studying the average. 
However, this opens the possibility of perturbing in e = 2 — d, as will be shown in this paper. 

The small-world network is constructed by starting with a regular lattice in c?-dimensions. 
Then, some set of long-range links are added: a given pair of sites i,j is connected with 
probability pa'^ /V, where V is the total number of sites in the system. Here, we define a 
length a as the lattice scale, and p as the density of links. Then, as — *■ oo, each site has 
a Poisson distribution of links emanating from it, with on average pa"^ links. Typically, the 
links, if any, leaving a given site will connect that site to other sites far away in the system. 

Looking for universal results, we consider the case of a low density of links, pa"^ << 1. 
Ignoring sample-to-sample fluctuations, the natural mean-field system to consider is one in 
which each site is coupled to all others with a strength ~ p/V. This leads to a solvable 
problem with a correlation length ^ (x p~^^^, or a correlation volume ^'^ oc p~'^^^. Returning 
to the original problem with fixed links, we see that such a volume has p^~'^/'^ links in it, and 
as p ^ 0, this number of links tends to zero for p < 2. This is a major problem. Mean-field 
theory ignores fluctuations in the number of links, which is only justifled if the number of 
links is large, a condition which is not satisfled in this case for d < 2. This problem is a result 
of a violation of the modified Harris criterion introduced in [5] . 

Instead, we expect that the correlation length for the average Green's function must be at 
least p~^^'^ for small p, so that there is on average at least one link in a correlation volume. We 
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will see below that the correlation length is in fact proportional to this limit, and the 

e-expansion will enable us to calculate the prefactor, as well as to study fluctuations about the 
average. Given a p~^/'^ scaling of the correlation length, the number of links in a correlation 
volume is some p-independent number. The basis for the e-expansion is the observation that, 
within a self-consistent approximation, this number diverges as 1/e, so that fluctuations about 
a self-consistent mean-field can be calculated diagrammatically as done below. 

Green's Function — Our general problem is to study diffusion on the small-world network. 
Thus, we must calculate Green's functions of the following problem: dtPi(t) ~ ~^j^ijPjj 
where Pi(t) is the probability of finding some randomly walking particle at site i and f is the 
Laplacian on the small-world network. We take 

r = r° + qu. (1) 

Here, F^^ is the Laplacian on the regular part of the network. = — a^^ if i and j are 
neighboring sites on the regular network, while T% is equal to a^^ times the coordination 
number of site i. Uij is the Laplacian on the long-range links. Uij = —1 if i,j are connected 
by a long-range link, where Uu is equal to the number of long-range links leaving site i. 

We have inserted an extra factor of q multiplying the matrix U in Eq. (Q. If g is small 
and p is large, this implies that we have a high density of weak links and the problem can be 
solved via a mean- field theory in which we ignore the fluctuations in the local density of links. 
However, we will be interested in the opposite case, where q is of order 1, while p is small. 

Fourier transforming, we are interested in the Green's function (iu! + F)"^. However, we 
will focus on the case w — > 0. The matrix F has a zero single mode, due to the conservation of 
J2i Pi by the diffusion process. Throughout, we will work in the subspace orthogonal to this 
zero mode, defining G — limi^^o(i^ + F)~^ and = lim^^o{iuj + F°)~^ in this subspace. 
This Green's function is related to the return probability of a random walker, and also to 
the roughness of a surface defined by Edwards- Wilkinson dynamics on the network [4]. We 
will compute Gy, where the line denotes averaging over the ensemble of different random 
networks, as well as computing higher moments such as GijGki- 

For p small (compared to a~'^), the length much larger than the lattice scale, so 

that we can take a continuum limit, setting F'' equal to the continuum Laplacian d^. In the 
limit of small p, the probability of a single site having more than one link becomes vanishingly 
small. This continuum limit leads to ultraviolet divergences for d > 2, which are cutoff at the 
lattice scale a. However, for d < 2, this limit is completely convergent, and thus for d < 2 
we find universal results, independent of the lattice details. This use of a continuum limit is 
essential to continue the results to arbitrary real dimension d < 2. 

For clearer notation, in the continuum limit we will often use d-dimensional vectors x, y, ... 
to label lattice sites, rather than indices .... We write G{x,y) to denote a Green's function 
between points x,y and G{x) to denote a Green's function G{x,0) (the averaging restores 
translational symmetry). We also use a vector k to label momenta, defining the averaged 
Green's function at Fourier mode k to be G{k) — j d'^xG{x) cxp(jfc • x). We use a similar 
notation for other matrices in the continuum limit. Discrete (5- functions get replaced by Dirac 
(5- functions, while sums get replaced by integrals. 

Ignoring fluctuations in the local density of links, the mean- field approximation consists 
of replacing U by U, so in the continuum limit G « {d^ + U)~^ . Then, we find that U{x^ y) = 
pqS{x — y) — pq/V. This gives a Green's function G(a;) = / d'^k{2TT)'~'^ exp{ik ■ x){k'^ + pq)~^ . 
The correlation length is {pq)~^^'^ as discussed above, and thus this expansion must breakdown 
for d < 2 in the limit of small p at fixed q. 

Self- consistent Calculation — To go beyond the mean-field calculation, we use impurity- 
averaged perturbation theory [6], following [4]. We continue with lattice notation here, since 
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a) = = - -^L^ 

b) E = X + (5^ + >^ + + 



Fig. 1 - a) Diagrammatics for E and b) Self-consistent sum of interactions with a single link. 

this development is not specific to the continuum limit. The perturbative expansion of G in 
powers of ?7 is G = — G^UG^ -t- G^U G^U — where a product of matrices is implied. 
Each of the terms in this expansion can be computed using the given distribution of disorder. 
The mean- field calculation is based on the following approximation: G « (F'^ + U)^^ = 
qo _ qOjjqo _|_ Q^uq^uq^ _ ... These two expansions are equal at zeroth and first order in 
J7, but differ at second order in U , by an amount depending on the second cumulant of U . 

It is useful to introduce a diagrammatic notation for the perturbative calculations, as 
shown in Fig. 1. We use a single solid line to denote G°, and a double solid line to denote 
G. We use a cross with no dashed hues attached to denote the average Uu. A pair of crosses 
connected by a dashed line is used to denote an average UuUjj — Uu Ujj. Three or more 
crosses connected by dashed lines are used to denote the third and higher cumulants. These 
diagrams denote two, three, or more scatterings of a single link. Circles are used to denote 
averages of off-diagonal terms, Uij for i j- Dashed lines can also connected both circles 
and crosses, again denoting higher cumulants. Within this diagrammatic approximation, the 
weight of a given diagram is equal to g)"^. Here n; is equal to the number of sets of 

circles or crosses connected by dashed lines, while is equal to the total number of circles 
plus crosses. Thus, each link gives one factor of p (the probability of finding a link) multiplied 
by —q to the number of times that link appears in U. 

We note that a single circle, not connected by a dashed line, denotes Uij, which vanishes 
as — > oo, and thus may be ignored in this perturbation expansion. In general, any diagram 
involving an odd number of off-diagonal terms for a given link vanishes as 1^ — » oo. 

Let us return to the continuum limit for specific calculations. We introduce the self-energy 
E(fc) by setting G(fc) = [F'^-|-E(fc)]~^, as shown diagrammatically in Fig. 1(a). Approximating 
E to leading order in U we get T,{k) = U(k). We have E(0) = 0, However, for all fc 7^ 0, E(fc) is 
/c-independent, so that, outside of the subspace of the zero mode, we may write E(fc) — T, — pq. 
We now compute corrections to this result, to second order in U. To this order, E(fc) is still 
fc-independent for k ^ 0. This self-consistent Born approximation is defined by the first three 
diagrams on the right-hand side of Fig. 1(b). Higher diagrams of Fig. 1(b) include all diagrams 
with arbitrary numbers of scatterings off a single link. We self-consistently use double lines 
for the Green's function in Fig. 1(b), giving T, = pq — 2pq^G{Q), where 




(2) 



For d < 2, the integral of Eq. (|2Jl is convergent, and equal to 



TT' 



.d/2 



(3) 



The self-consistent Born approximation at second order in U then becomes E — pq — 2pq^G{0). 
As a first guess at solving this equation, we substitute the first order result for E in the equation 
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for G(0), getting T, ^ pq - 2p'^/'^q'^/^+^TT'^/^T{l - d/2)/{2TTY for d < 2. The second order 
correction is comparable to the first when q ^ p2/d-i q;"i^^s^ fQj- q <<; pZ/d-i^ perturbation 
theory may be appHed. We wiU be interested in the opposite hmit of p small at fixed q so 
that q » p'^l'^~^ and perturbation theory breaks down for d < 2. 

A better approximation for S is to include all diagrams in Fig. 1(b), summing all diagrams 
involving interactions with a single link. Since we are interested in the case in which the density 
of links is small, but the scattering q off a single link is of order unity, we sum all scatterings 
off a single link. This gives E = pq/[l + 2qG{0)]. For p small, we will find G{Q) >> 1, so that 
we can take E — pq/[2qG{0)] = p/[2G{Q)]. Here, the factor of two arises from the different 
possible diagonal and off-diagonal scatterings off a single link: if a particle interacts with a 
given link n times, then there are 2"~^ total diagrams contributing to S. Using Eq. for 
G(0), we find for d < 2 that the solution of E = p/[2G(0)] is E = 

Eo = [(2^) V(27r''^'r(f - d/2)]2/''. (4) 

We have placed the subscript on E, since later this result will be used as a zero order ap- 
proximation, with corrections to it in powers of e. This gives E oc p^/"*, so that the correlation 
length varies as p^^"^. We will find below that a more careful treatment leads to universal 
corrections in orders of e to the prefactor of Eq. but does not change the scaling with p. 

In contrast, for d > 2, the integral of Eq. ||2Jl is divergent at large k, but it converges at small 
k even for E = 0. The divergence at large k is cut off by the lattice scale, a. Thus, G(0) has 
a well-defined limit as E — > 0; this limit has a non-universal dependence on the lattice details 
and is equal to some number, g. The self-consistent equation becomes — pq/[l + 2qG{0)]. 
To leading order in p, the solution of this equation gives E = pq/{l + 2qg). The magnitude 
of the corrections to the mean-field result, E — pq, depends on the product qg and does not 
have any universal behavior. There are also further corrections to the mean- field result which 
involve scattering off multiple links. These corrections will not be considered here. 

Finally, consider d = 2. In this case, we will show below that the approximation above 
leads to exact results for the small p behavior of E. The integral of Eq. ((2Jl is logarithmically 
divergent. The divergence is cut off at a A: of order a~^, giving G(0) = — (27r)~^ log(aE-^/^). 
The self-consistent equation then gives E = pq/[l + 2gG(0)]. As for d < 2, G(0) >> I for 
small p so that this reduces to E = p/[2G(0)] = — 27rp/ log(a^E). Solving this self-consistent 
equation for a'^p « 1 gives E = — 27rp/ log(pa^), plus subleading terms. 

Corrections in Two Dimensions and Renormalization Group — Thus far, the calculation 
has followed [4]. Now, however, we go beyond the self-consistent calculation to obtain an 
expansion in e. The first step is to analyze the behavior for d = 2 in more detail, and show 
that the results above are exact for the leading scaling of E with p. In the next section, a 
diagrammatic perturbation expansion will be presented to compute results in powers of e. 

In d = 2, the correlation volume is of order 1/S. The physical motivation for the results 
below is that for E ^ —p/ log(pa^), the average number of links in a correlation volume, p/E, 
is of order — log(pa^). Thus, the average number of links in a correlation volume diverges as 
p — > 0. This, however, is exactly the condition we need to make mean-field theory work, as it 
enables us to ignore fluctuations in the number of links. 

Turning away from perturbation theory, we now consider the problem scale-by-scale, start- 
ing with the shortest distances. We begin with a physical description of the problem, followed 
by a more careful calculation below. On length scales I << p~^/'^, each link can be considered 
independently, since the density of links at such scales is very small: the probability of finding 
a link near any other link is negligible. Solving the problem of a single link at a scale I is 
very similar to the calculation done above. We must sum multiple scatterings off the single 
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)(=)('* Q^' + )f^>c'+ 
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X 

Fig. 2 - a) Definition of resummed scattering, b) Disallowed diagrams, c) Leading corrections to G. 
d) Leading contribution to GG — GG. For all diagrams, single line denotes G' in this figure. 

link. This will renormalize the interaction strength of the link from q to q/[l + 2qGi{Q)], 
where G;(0) ~ (27r)~^ log(^/a) is the Green's function of matrix r° at scale I. Any function 
which cuts off momenta less than l^^ will suffice for G;. We choose Gi{k) — (fc^ + . 
Then, once we reach scale / = the link interaction strength has been renormalized to 

qren = q/[l + 2qG,(0)] - -27r/log(pa2). 

The number of links in an area of size a? was pa^ . In contrast, on the new scale I = p^^/"^ , 
the number of links in area P is of order unity, while the scattering strength qren is small, 
and thus corrections to mean-field theory become negligible for small p. Then, we find S = 

Pqren = -2'Kp / \og{pa'^) . 

Expansion in e — Now, we wish to compute quantities for d < 2. In this case, we find 
that once we reach a scale I ~ p~^/'^, the link interaction strength qren has been reduced to 
an amount of order q/[\ ~\- 2gG/(0)] ~ el^'^. Unlike the case in two dimensions, qren does not 
vanish as we take the limit pa'* 0. However, we will find that the interaction strength of 
the links is of order e, so that corrections to mean-field theory can be computed in powers of 
e. The e-expansion is defined by a simple re-ordering of diagrams. First, note that Eq. Q 
defines a length scale I by 

/ = £(^1/2 ^ [(27r) V(27r'*^'r(l - d/2)Y^/'^. (5) 

so Gi{k) — (fc^ -I- Then, we define the resummed diagonal scattering off a single link 

as in Fig. 2(a), denoted by a cross with double dashed lines. Here, the single line denotes 
the Green's function G;, rather than G° as before (this resummation is closely related to the 
T-matrix in scattering theory). The resummed off-diagonal scattering is defined similarly by 
a circle with double dashed lines. 

The starting point for the e-expansion is the self-consistent calculation of Eq. Q. That 
expansion resums a large set of diagrams known as the "rainbow" diagrams; this is the set 
of all diagrams except those in which a particle first scatters off one link, then off a second, 
then returns to the first link, and then returns to the second link (possibly with additional 
scatterings off other links included) . We now present an expansion for G which exactly includes 
all of the missing diagrams, such as in Fig. 2(c), with the correct coefficent, using resummed 
scatterings with each link. We then show that this leads to an expansion in powers of e. 

First, write down all possible diagrams, using G/ for the Green's functions and double 
dashed lines for interactions with links, subject to the following two constraints: (1) dou- 
ble dashed lines must always connect at least two crosses or circles. That is, single crosses 
or circles never appear individually, so the first diagram of Fig. 2(b) is not allowed. This 
constraint on diagrams is included because such diagrams are already taken into account in 
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the self-consistent expansion, by the replacement of G° by G;. (2) two successive interactions 
with links must always involve different links (links not connected by dashed lines), so that the 
second diagram of Fig. 2(b) is also not allowed. The second constraint is presented to avoid 
overcounting multiple scatterings. Given these constraints, the leading correction to the aver- 
aged Green's function is given by Fig. 2(c). Each diagram is assigned a factor of p"' (— qren)"° , 
where ni is now the number of sets of circles and crosses connected by double dashed lines 
and Us is again the total number of circles and crosses. Here, qren = + 2<7G;(0)] = So/p. 

This expansion, though it simply reorders the diagrams, perfectly realizes the desired 
treatment of the problem scale-by-scale. We must start with the shortest scales. However, at 
the shortest scales, the only processes involve multiple scattering off single links; only when 
the scale becomes of order 'p'^^l'^ do processes with multiple links become important. Thus, 
we resum scattering off the single link, up to the scale I. Unlike the case for d = 2, processes 
involving multiple links, such as that shown in Fig. 2(c), are no longer vanishing, but we now 
show that they are higher order in e. First, it is convenient to rescale all distances by /. Then, 
we set X = x/l, k = kl. We use the Green's function Gi(fc) = (fc^ -I- 1)~^, and define p — pl"^, 

while Qren = Qrenl^''^- Then, wc define G(i) by G(x) = G{x)l'^~'^. Then, the perturbation 
expansion for G is obtained by using the same set of diagrams as above, but replacing p by p, 
Qren by qren, and Gi by Gi. Now, we have 

p = 27r'*/2r(l - d/2)/(27r)'^ = e-HTT'^^y{2TTf + (6) 

where here we have given a series expansion of the results for p, qren, and the ... denote higher 
order terms in e. Thus, we have defined a new problem of scattering with links of strength e 
and density e~^; this a problem of a high density of weak links so that perturbative techniques 
work well and lead to an expansion in e. Physically, one can imagine that for d close to two, a 
randomly diffusing particle has only a weak interaction with a link: the dimension of the path 
of the particle is two, so that for d close to two the particle can easily "miss" a given link. 

There are a finite number of diagrams at each order in e as we now show. The number of 
scatterings of each links is at least two by the rules above, so > 2n/. Then, p"'g"g„ ^ ^n^-ni 
is at least order e"'. Thus, to order e", we need only consider diagrams ni < n. For a given 
ni , we need only consider diagrams with Ug < n + ni, leaving us with only a finite number of 
diagrams at each order in e. It is important to consider the possibility of ultraviolet divergences 
in two dimensions, since they may give extra factors of e~^. However, the only ultraviolet 
divergences in two dimensions in this expansion arise from Green's functions Gi (0) and this 
expansion resums all such divergent diagrams which involve scattering off at least one link: 
the only divergent contribution to G(0) is from the very first diagram in Fig. 2(c). Finally, to 
all orders in e, the scaling of G(0) with p is unchanged from the self-consistent calculation. 

We now use this formalism to compute specific results. We first compute the average 
Green's function. From Fig. 2(c), we have G(0) = p-'^[Gi{0)+p^q^^„ J J d'^xd'^yGi{x)Gi{y)Gi{x- 
uY + ■•■]■ We have numerically evaluated this integral in rf = 2 to get the leading corrections 
in e (to simphfy this integral, we used the trick that it equals — (1/4)9; / d^a:G/(x)^ at I = 1). 
The resuh is G(0) = p-'^[Gi{G) + e27r2(.00106...) + ...], where we have also evaluated Qren 
for d = 2, and where I is given by Eq. ©. Here, Gi(0) tt'^/'^ / {2ttYT{1 - d/2) and hence 
diverges as e~^. As noted above, this is the only ultraviolet divergent diagram. 

To go to higher order in e, it would be necessary to keep additional diagrams, as well as 
to expand the integrals in powers of e near d =2. Since we are interested in d = 1, let us first 
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evaluate Fig. 2(c) in c? = 1, where p = q^en = 1, so 

G(0) =P"MGi(0) + j J dxdyGi{x)Gi{y)Gi{x-yf + ...]= p-\l/2 + 5/256 + ...). (7) 

At higher order, if we evaluate all diagrams in o? = 1, but use the e expansion to define the 
ordering of diagrams for us, all the diagrammatic integrals can be performed exactly, since 
Gi{x) then has a simple exponential decay. This is a task left for future work, as is a test of 
the convergence properties of this expansion. Eq. {Tj) compares well with numerics [4]. There, 
G(0) was found to scale p~^, with a prefactor slightly larger than 1/2, close to the result here; 
comparison of the exact difference between the prefactor and 1/2 would require higher order 
calculations here and larger system sizes in the numerics. 

Fluctuations in the Green's function can also be computed using these techniques. From 
Fig. 2(d), we have G(0)2 - (G(0))^ = P~^U dxGi{x)'^ + ...) = ^"2(1/32 + ...) 

Discussion — We have presented an e-expansion for the properties of the small-world net- 
work in d = 2 — e dimensions, enabling us to compute averages of moments of the Green's 
function. This technique should, however, have much greater generality. In [5], general criteria 
were put forth for the breakdown of mean-field theory for any statistical mechanical system 
on a small- world network, in analogy to the usual Harris criterion for disordered systems [7]. 
Hopefully, for other systems in which these criteria are violated, it will be possible to pro- 
vide an e expansion near the critical dimensionality, analogous to what has been done for the 
violation of the ordinary Harris criterion [8] . 

These criteria can be extended to other problems than networks, such as a randomly 
diffusing particle in an array of traps [9], a model inspired by work in reaction-diffusion 
processes [10]. The criteria [5] correctly predict d = 2 as the region for the breakdown of 
mean-field theory, and the techniques here can be used to compute average properties of the 
Green's function in d = 2. For the random trap system in d = 1 with p << q, the particle has 
only a small probability to pass through one trap to reach the next trap which is typically far 
away. This makes the trap system very far from mean-field theory and suggests that the e- 
expansion will not converge to e = 1; however, it also enables the solution of the d = 1 system 
by considering each interval between traps independently. In contrast, for the network, the 
particle can interact with several links, even for p —> 0, thus suggesting that the perturbative 
e-expansion may work for the small-world case. A true test will require evaluation of additional 
higher order diagrams and is left for future work. 
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